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Abstract 

We solve N = (8,8) super Yang-Mills theory in 1+1 dimensions at strong 
coupling to directly confirm the predictions of supergravity at weak coupling. We 
do our calculations in the large- N c approximation using Supersymmetric Discrete 
Light-Cone Quantization with up to 3 x 10 12 basis states. We calculate the 
stress-energy correlator (T ++ (r)T ++ (0)) as a function of the separation r and 
find that at intermediate values of r the correlator behaves as r -5 to within 
errors as predicted by weak-coupling supergravity. We also present an extension 
to significantly higher resolution of our earlier results for the same correlator in 
the M = (2, 2) theory and see that in this theory the correlator has very different 
behavior at intermediate values of r. 



1 Introduction 



We give a direct evidence of the duality between supersymmetric field theory at strong 
coupling and supergravity at weak coupling. This conjectured duality, known as the 
Maldacena conjecture [1], has been studied extensively, and there is considerable qual- 
itative support for it. Most, if not all, of these tests are in some way indirect, because 
of the difficulty of solving supersymmetric field theories at strong coupling. The result 
we present here is unique, in that it is based on a numerical solution of an exactly 
supersymmetric field theory at strong coupling. 

Our method is Supersymmetric Discrete Light-Cone Quantization (SDLCQ) [2, 3]. 
It is a well established tool for calculations of physical quantities in supersymmetric 
gauge theory and has been exploited for many supersymmetric Yang-Mills (SYM) theo- 
ries. Briefly, the SDLCQ method rests on the ability to produce an exact representation 
of the superalgebra but is otherwise very similar to Discrete Light-Cone Quantization 
(DLCQ) [4, 5]. In DLCQ, we compactify the x~ = (t — z)j\p2 direction by putting 
the system on a circle with a period of 2L, which discretizes the longitudinal momen- 
tum as p + = nn/L, with n a positive integer. The total longitudinal momentum P + 
becomes Kn/L, where K is an integer known as the harmonic resolution [4]. The 
positivity of the light-cone longitudinal momenta then limits the number of possible 
Fock states for a given K, and, thus, the dimension of Fock space becomes finite, en- 
abling us to do some numerical computations. It is assumed that as K approaches 
infinity, the solutions to this large, finite problem approach the solutions of the field 
theory. The difference between DLCQ and SDLCQ lies in the choice of discretizing 
either P or Q~ to construct the matrix approximation to the eigenvalue problem 
M 2 \m) = 2P+p-\^} = 2P+(Q-)7v / 2~|^), with P+ = Kn/L. For more details and 
additional discussion of SDLCQ, we refer the reader to Ref. [3]. 

The M = (8,8) supersymmetric theory in 1+1 dimensions in the large- N c limit is 
discussed in Refs. [6, 7, 8]; however, the published results are primitive compared to 
what can be obtained today, because of our greatly improved hardware and software. 
In the present paper we are able to present results for a resolution of K — 11. At this 
resolution the problem has 3 x 10 12 basis states. By fully exploiting the symmetries of 
this theory, we only need to deal with a subset of 3 x 10 7 such states. 

We will look at the two-point correlation function of the stress-energy tensor, namely 
(T ++ (r)T ++ (0)). The expected behavior in the ultraviolet (UV) and infrared (IR) 
regions is 1/r 4 , and the predicted behavior in weak-coupling supergravity theory in the 
intermediate region is [9] 1/r 5 . We find the power behavior in M = (8, 8) theory in the 
intermediate-r region to be consistent with 1/r 5 . 

An interesting effect in the calculation is that the finite-dimensional representations 
in SDLCQ with odd and even values of K result in very distinct solutions of the SYM 
theory. We saw this result in our work on M = (2, 2) SYM theory in two dimensions [10] . 
Only at infinite resolution are the solutions identical. One might initially think that 
having two numerical representations of the supersymmetry is a shortcoming of the 
SDLCQ approach, but it turns out to be an advantage because it provides an internal 
measure of convergence. We present new results for the M = (2,2) theory, up to 
resolution K = 14, which clearly demonstrate this. 

These results for N = (2,2) also correct a sign error in our earlier numerical results 
for Af = (2, 2). The quantitative change is very small and has no qualitative effect. 
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The structure of this paper is the following. In Sec. 2 we review A/"=(8,8) SYM 
theory. In Sec. 3 we discuss the symmetries of this theory, which play a critical role in 
allowing us to greatly reduce the basis set we need to retain. We discuss the two-point 
correlation function of the stress-energy tensor in Sec. 4. Finally, in Sec. 5, we present 
our numerical results. A summary and some additional discussion are given in Sec. 6. 



2 Review of M = (8, 8) SYM Theory 

Before giving the numerical results, let us quickly review some analytical work on 
Af = (8, 8) SYM theory, for the sake of completeness. For more details, see Ref. [6]. This 
theory is obtained by dimensionally reducing Af = 1 SYM theory from ten dimensions 
to two dimensions. In light-cone gauge, where A_ = 0, we find for the action 



= J dx + dx~tr 



d+Xjd-Xj + ie T R d + e R + ieTd-e L 



,2 



+\{d_A + Y + gA + J+ + V2g0 T MXj,e R ] + ^[Xj,Xj] 2 

where x ± are the light-cone coordinates in two dimensions, the trace is taken over 
the color indices, the Xj with / = 1, . . . , 8 are the scalar remnants of the transverse 
components of the ten-dimensional gauge field A^, the two-component spinor fields 9 R 
and 6l are remnants of the right-moving and left-moving projections of the sixteen- 
component spinor in the ten-dimensional theory, and g is the coupling constant. We 
also define the current J + = i[X/,<9_Y/] + 29 R 9 R and use the matrices $ given in 
Ref. [11]. 

After using the equations of motion to eliminate all the non-dynamical fields, we 
find for P~ = J dx~T + ~ the expression 

P- = 9 2 J dx-tr {-^J + ^J + ~ \[Xi,Xj? + ^[XiMF^PAXjJr]) ■ (2) 

The supercharges are found by dimensionally reducing the supercurrent in the ten- 
dimensional theory. We find 



Qa=gJ dx-tr (^-2^J + ^-u a + 2- 1 /H[X I ,X J ]((3 I (3j) ai] u^ , 



(3) 



where a, r] — 1, . . . , 8 and the u a are the components of 9 R . We expand the dynamical 
fields Xj and u a in Fourier modes as 

Xi Pq (x~) = -= J -^=[A Ipq (k + )e- lk+x - + A\ qp (k+)e ik+ n, (4) 



2tt ^2k A 

j roo rth^ , , 

u apq (x-) = ^J o -^[B apq (k+)e^ + *- + Bl qp (k + )e ik+ *-}, (5) 

where p,q = 1,2, . . . ,N C stand for the color indices, and A and B satisfy the usual 
commutation relations 

[A Ipq (k + ),4 rs (k' + ]] = 5 IJ 5 pr 5 qs 5(k + -k ,+ ), (6) 

{B apq (k + ),Bl rs (k' + )} = 5 a(3 5 pr 5 qs 5(k + -k' + ). (7) 
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We work in a compactified x~ direction of length 2L and ignore zero modes. With 
periodic boundary conditions, momenta are restricted to a discrete set of values [2] 
k + = irk/L, with k a positive integer. The integrals over k + are replaced by sums: 
/ dk + — > j SfcLi, and Dirac delta functions become Kronecker deltas: 5(k + — k' + ) — > 
^8 kk t. We then rescale the annihilation operators 

^a(k) = A(k+ = ^), ^b(k) = B(k+ = ^), (8) 

so that 

[a Ipq (k),a]j rs (k')] = 5 u 5 pr 5 qs 5 kk >, {b apq (k), b^ rs (k')} = 5 a/3 5 pr 5 qs 5 kk >. (9) 
The expansions become 

Xi Pq (x~) = -j= J2—[a Ipq (k)e- i i kx ~ + a\ qp (k + )e^-}, (10) 

1 00 1 

u am {x-) = —^—[b^We^ + b^Wji**-]. (11) 



In terms of a and b, the supercharge is given by 



oo 



Qa = — \/§ Yl 6 (ki+k2),k 3 \ ( 12 ) 



' ^ /Ci,fc2,fc3 = l 



X 



^ [^Lj( /c 3)a/i m (A;i)a /mi (A; 2 ) - a} im (fci)4 mj - (^2)^(^3)] 



+ ( fel ^ fe3 ) [ fl ^( A; l) & Lj(^2)a/ii(A;3) -a/ij( fc 3)a/im(fcl)6amj(fc2)] 

+ VK (~^T^) ^(^^(^W^) - &Lm(fci)aLj( fc 2)aj«(A;3)] 

-^-K^(fe)fcmm(A:i)6^(A;2) + ^m^O 6 ^* 2 ) 6 ™^)] 

_ r ( fc 3)6vm(A;i)6amj (fo) + &Ji m (*l)&Lj (h)b v ij (fc 3 )] 
#2 

+7- (h)b V i m (ki)br,mj (k 2 ) + 6j im (A;i)6j roj -(A;2)6aij(A;3)] 
+2^ 4 ^-^- [&Lj(^3)aj»m(fei)ajmi(fc2) + a/ im (A;i)a} mj (A;2)6 ai j(A;3)] 
+ 4 ^g^ [ofij(fe3)&aim(fei)a fm j(A;2) + 6Lm(^i)aLj( fc 2)a/ij(A;3)] 

+ ^= [Q|ij(^3)Qjim(fel)&amj(fe2) + 4im(fcl)&Lj (faKij (*%)] j |, 

where we have used the relation ([/?/, flj]) a ,ri = 8 ari 5ij. 
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3 Symmetries 



The superalgebra relations that involve Q+, as specified by the anticommutators 

{Q+,Q+} = 5 a/3 2 v / 2> + , {Q+,Q^} = 0, (13) 
are satisfied, but the anticommutators for Q~ 

{Q-,Qp} = 5 a p2V2P~ (14) 

become in SDLCQ 

{Q-,Qp} ^0 if (3, {Q~,Q~} = 2V2P-. (15) 

Although we have different P~ for different Q~, we can find unitary, self- adjoint trans- 
formations C a ' a , such that C a i a P~C a i a = P~,. Thus the eigenvalues of the different 
P~ are the same, and we may choose any one of the Q~'s, at least for our purposes. 
In what follows we will use Q$ and will suppress the subscript unless it is needed for 
clarity. 

To reduce the size of the numerical calculation, we seek symmetries that block 
diagonalize the P 8 ~ matrix. One such symmetry is a Z2 symmetry of Q 8 , called S- 
symmetry [12]: anj — > —ajji, b a ij — > —b a ji. Extending the work of [8], we also look for 
permutations and sign changes of the bosonic and fermionic operators that leave Q 8 
unchanged. The terms not involving (3 matrices are preserved under all interchanges 
that fix b 8 . The additional condition is that (PiPj — PjPf) a /3b ji i 3 ajaj be invariant with 
respect to the replacement 

ai -> r(/)ap(/), bp -> s{f3)b q{f}) . (16) 

Defining N a (f3, 1, J) = (PiPj — PjPf) a /3, this gives us the condition 

N a (P, I, J)b\a ia j = s(p)r(I)r(J)N a (P, J, J)b\ w a p{I) a p{J) . (17) 

Each possible flavor combination appears at most once on each side, so that we can set 
the coefficients equal, term by term, and obtain the condition 

N a (q([3),p(I),p(J)) = s([3)r(I)r(J)N a (P,I,J). (18) 

This is then used to determine the allowed transformations as specified by r, s, p, and 
q. There are (8!) 2 2 16 such transformations. When we restrict these to transformations 
that leave b 8 unchanged, and thus have q(8) = 8 and s(8) = 1, we have instead a set of 
8!7!2 15 w 6.7 x 10 12 transformations. 

To determine the allowed transformations, we first find the permutations such that 

\N a (q(P),p(I),p(J))\ = \N a (P,I,J)\. (19) 

Then, among these we check for choices of r(J) and s(I) that allow the full condition (18) 
to be satisfied. There are 1344 permutations that satisfy the absolute- value condition 
(19), and, for each of these, there are 16 choices of r(J) and s(I) that satisfy the full 
condition (18). From these, we find that the group of transformations that leave Q$ 
unchanged can be generated by 7 operators that square to 1. These generators are 
listed in Table 1. 
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Table 1: Generators of transformations that leave Q$ unchanged. The first three 
transformations are permutations with sign changes, and the last four involve only sign 
changes. 
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be 
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b 2 
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ai 
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-67 


4 


-Ol 


-ai 
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b 5 
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ai 


a 2 


a 3 


— 04 


-05 


a 6 


-a 7 


-a 8 


61 


fo 2 


-63 
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-be 


6 7 
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a 2 


03 
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a 5 


-a 6 


-a 7 


a 8 
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-62 


63 


-64 


-b 5 


be 


-67 
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ai 


-a 2 


a 3 


a 4 


-a 5 


-o 6 


-a 7 


a 8 




fo 2 


63 


-bi 


b 5 


-be 


-67 



4 Correlation functions 

One of the physical quantities that we can calculate nonperturbatively is the two-point 
function of the stress-energy tensor. Previous calculations of this correlator in this and 
other theories can be found in [7, 8, 13]. Ref. [8] gives results for the theory considered 
here but only for resolutions if up to 6. We can now reach K — 11, with its 3 x 10 12 
basis states, by careful use of the symmetries discussed in the previous section. The 
largest matrices we have to consider are 3 x 10 7 by 3 x 10 7 . 

We find that there are distinct behaviors in the correlation function for even and 
odd K. This is a familiar aspect of SYM theories with extended supersymmetry, and 
we have argued [10] that we have two different classes of representations at finite K, 
which become identical as K — > 00. 

Let us first recall that there is a duality that relates the results for the two-point 
function in M = (8,8) SYM theory to the results in string theory [8]. The correlation 
function on the string-theory side, which can be calculated with use of the supergravity 
approximation, was presented in [7], and we will only quote the result here. The 
computation is essentially a generalization of that given in [14, 15]. The main conclusion 
on the supergravity side was reported in [9]. Up to a numerical coefficient of order one, 
which we have suppressed, it was found that 

AT 3 / 2 

(T++(r)T++(0)> = — £-^. (20) 
Qym^ 

This result passes the following important consistency test. The SYM theory in two 
dimensions with 16 supercharges has conformal fixed points in both the UV and the IR 
regions, with central charges of order and N c , respectively. Therefore, we expect the 
two-point function of the stress-energy tensor to scale like iV c 2 /r 4 and N c /r 4 in the deep 
UV and IR regions, respectively. According to the analysis of [16], we expect to deviate 
from these conformal behaviors and cross over to a regime where the supergravity 
calculation can be trusted. The crossover occurs at r = 1/ QymVNc an d r = ^/W c / 9ym- 
At these points, the N c scaling of (20) and the conformal result match in the sense 
of the correspondence principle [17]. We should note here that this property for the 
correlation functions is expected only for M = (8, 8) SYM theory. However, it would 
be natural to expect some similarity between Af = (8, 8) and Af = (2, 2) theories. 
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We wish to compute the expression F(x~, x + ) = (T ++ (x~ , x + )T ++ (0, 0)), where we 
fix the total momentum in the x~ direction. It is more natural to compute the Fourier 
transform and express the transform in a spectrally decomposed form [7, 8] 

P(P_,:r+) = -L(T++(P_,:r + )T++(-P_,0)> (21) 

= E^(0|T ++ (P_,0)|.)e-^ + (.|T ++ (-P_,0)|0). 

The position-space form of the correlation function is recovered by Fourier transforming 
with respect to P_ = P + = Kir/L. We can continue to Euclidean space by taking 
r = \/2x + x~ to be real. The result for the correlator of the stress-energy tensor was 
presented in [7], and we only quote the result here: 



7T 




2 



M: 



4 



F(x-,x+) = Y,t(0\T ++ (KM (-) -J^K^y/to^F), (22) 



where Mj is a mass eigenvalue and K±(M(r) is the modified Bessel function of order 4. 
In [7] we found that the stress-energy operator is given by 



T ++ (x-,x + ) = tr 



(d_X 1 ) 2 + Uiu a d_u a - i(d_u a )u a ) 
2 



(23) 



When written in terms of the discretized creation operators, at and b\ as in Eqs. (10) 
and (11), we find 



7T 



T ++ (K)\0) = ^E[-^k(K-k)a\ tl (K-k) a yk) (24) 



+ (f - k) bl :j (K - k)bi n (k) 



|0). 



The matrix element (L/7r)(0\T ++ (K)\i) is independent of L and can be substituted 
directly to give an explicit expression for the two-point function. We see immediately 
that the correlator behaves like 1/r 4 at small r, for in that limit, it asymptotes to 

On the other hand, the contribution to the correlator from strictly massless states is 
given by 




L 



7T 



2 6 



Eb^MW^Sw- (26) 



That is to say, we would expect the correlator to behave like 1/r 4 at both small and 
large r, assuming massless states have non-zero matrix elements. 

The operator used for calculating the correlator, T ++ (K) as given in Eq. (24), 
is preserved under all of the Z 2 symmetries discussed in Sec. 3. Therefore, we only 
need to consider states in the singlet sector. Including factors of 2 for supersymmetry 
and the S symmetry, the dimension of the Q~ matrix is thus reduced by a factor of 
86016. Because the dimension of the basis increases by a factor of almost 17 for each 
unit increment in K, symmetry allows us to increase the resolution by more than four 
without increasing the size of the matrix to be diagonalized. 
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5 Numerical results 



To compute the correlator using Eq. (22), we approximate the sum over eigenstates 
by a Lanczos [18] iteration technique, as described in [8, 13]. The results are shown 
in Figs. 1 and 2, which show the log-log derivative dlog 10 (/)/dlog 10 (r) of the scaled 
correlation function 

^(^-.^(O))^ 2 ^^. (27) 




Figure 1: The log-log derivative dlog 10 (/)/dlog 10 (r) of the scaled correlation function 
/, as defined in Eq. (27) of the text, versus log 10 (r), for the J\f = (2,2) theory. The 

separation r is measured in units of ^Jn/g 2 N c . In (a), the lines correspond to different 
values of resolution K from 3 to 14, with dashed lines for even K and solid for odd. 
The darker lines are for K = 13 and K = 14; the lower-.fr lines converge to these two. 
In (b), the vertical bars span four extrapolations to infinite resolution obtained from 
separate quadratic and cubic fits to even and odd K. 

Let us discuss the behavior of the correlator at small, large, and intermediate r. 
First, at small r, the graphs of / for different K approach as K increases. This 
follows Eq. (25), which gives the form / = 1 — 1/K. Second, at large r, obviously, the 
behavior is different for even and odd K . However, the difference gets smaller as K 
gets larger. The reason for this difference is that there are different finite-dimensional 
representations of the super-algebra [10] for even and odd K, which only coincide 
at infinite resolution. Looking at the detailed information of the computation of the 
correlator at larger r, we found that for even K there are massless states that contribute 
to the correlator, while there is no massless state that makes any contribution for odd 
K. Instead, it is the lowest massive states that contribute the most for odd K, and 
these states become massless as the resolution approaches infinity. 
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0.5 

log r 



0.5 

log r 



(b) 



Figure 2: Same as Fig. 1, except for the Af = (8,8) theory with K = 10 and 11 the 
largest values of the resolution. 

The results for the Af = (2,2) theory are shown in Fig. 1, and the results for the 
Af = (8, 8) theory in Fig. 2. The results for the Af = (2, 2) theory are obtained for 
higher resolution and more clearly show the convergence of the two representations. We 
are now able to extend our calculation of the Af = (2, 2) SYM theory from K = 11 [10] 
to K = 14. To calculate continuum results in SDLCQ, it is customary to extrapolate 
the finite resolution results to infinite resolution. We have done this for the curves for 
the Af = (2, 2) and (8, 8) theories by performing the extrapolation at successive values 
of r. We fit the even and odd resolution points with separate quadratic and cubic 
functions of 1/K to perform four extrapolations. Typical results of these calculations 
are shown in Fig. 3 for the Af = (8,8) theory at log 10 (r) = 0.2 and at log 10 (r) = 0.5, 

where r is measured in units of ^Jn/g 2 N c . 

In the intermediate-r region for the Af = (2, 2) theory, clearly the behavior of the 
correlation in Fig. 1(a) changes, and there is a flat region. The behavior is even more 
apparent in Fig. 1(b), the extrapolation to infinite resolution. In [10] we first saw the 
flat region for the Af = (2, 2) SYM theory, indicating that the correlator in this case 
behaves like 1/r 4 ' 7 . Note that the region of flattening extends farther out as K gets 
larger, for both odd and even K, implying again that the even and odd representations 
appear to agree as K goes to infinity. Unfortunately, there is currently no known 
prediction for the behavior of the correlator in the intermediate-r region this theory 
or even any reason to believe that the correlator should have a special behavior at 
intermediate values of r. 

In the intermediate-r region for the A/"=(8,8) theory, there is a prediction that the 
correlation should have a special behavior. This predicted behavior is very different 
from the behavior we found in the Af = (2, 2) SYM theory. We expect from Eq. (20) 
that the behavior is 1/r 5 . In [8] we found that the correlator may be approaching 
this behavior, and we indicated that conclusive evidence would be a flat region in the 
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dlog(/) 




dlog(r) 
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1/K 1/K 
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Figure 3: Sample quadratic (solid) and cubic (dashed) fits for extrapolations to infinite 
resolution in the J\f = (8,8) theory. The log-log derivative d\og 1Q (f) / d\og w (r) of the 
scaled correlation function /, as defined in Eq. (27) of the text, is plotted versus 1/K 
for (a) log 10 (r) = 0.2 and (b) log 10 (r) = 0.5, where r is measured in units of ^Jir/g 2 N c . 
Computed values of the log-log derivative are marked by circles for even K and triangles 
for odd K. 

derivative of the scaled correlator at a value of —1. In Fig. 2(b) we see that as r 
increases, the slope approaches —1. It overshoots —1 and then appears to approach 
— 1 from below and to remain consistent with that value for an extended range of r. 
Beyond this point, we see that the errors increase significantly. Thus, there is a finite 
range of intermediate values of r over which the numerical solution of strong-coupling 
SYM theory is consistent with the 1/r 5 prediction of weak-coupling supergravity. 

6 Discussion 

We have presented numerical results for the two-point correlation function of the stress- 
energy tensor, using SDLCQ in 1 + 1 dimensions in the large-iV c approximation, for 
J\f = (2, 2) SYM theory up to resolution K = 14 and for J\f = (8, 8) SYM theory up 
to resolution K — 11. There are two distinct classes of representations for these SYM 
theories, one for the resolution K even and one for K odd, and these representations 
become identical as K — > oo. The two-point correlators of the stress-energy tensor 
behave like 1/r 4 in the UV (small r) and IR (large r, K even) regions. The large-r 
behavior for K odd, on the other hand, has an exponential decay, since in the odd K 
representations the massless states become massless only as K — > oo. 

While there is no known prediction for the correlator of the stress energy tensor at 
intermediate values of r for the Af = (2, 2) theory, it is interesting that we find the 
dominant power law behavior 1/r 4 ' 7 . There is clear convergence over a wide range of 
intermediate values of r. 
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The results presented here also include a correction to our earlier work on the 
M = (2, 2) theory. We recently found a sign error in one term in our numerical calcu- 
lations. This correction does not change the qualitative behavior found in our previous 
calculation; the quantitative change is quite small. 

In M = (8, 8) SYM theory in 1 + 1 dimensions, the correlator is expected to behave 
like 1/r 5 in the intermediate region. We were able to confirm this power-law behavior 
with a flat region in the derivative of the scaled correlator. We have done this by 
calculating the correlator at successive values of the resolution and extrapolating the 
function to infinite resolution. The maximum resolution we reached was K — 11. At 
that resolution the total number of basis states in the SDLCQ approximation is 3 x 10 12 . 
Using all of the symmetries of the M = (8,8) theory, including some that are specific 
to the sector where the correlator is non-vanishing, we were able to reduce the number 
of basis states that we needed to consider to 3 x 10 7 . 

While just a few years ago it seemed inconceivable that one could reach these reso- 
lutions in SDLCQ, advances in our understanding of the problem as well as advances 
in software and hardware have allowed us to do so. We remain confident that we can, 
in fact, reach even higher resolutions. At this point we believe that the most important 
direction is to increase the number of fields that we consider. In particular, we are 
interested in considering theories with both extended supersymmetry and fundamental 
matter. 
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